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Abstract 

Over the last decade, implementations of several desingularization al- 
gorithms have appeared in various contexts. These differ as widely in their 
methods and in their practical efficiency as they differ in the situations in 
which they may be applied. The purpose of this article is to give a brief 
overview over a selection of these approaches and the applicability of the 
respective implementations in Singular. 

1 Introduction 

The problem of resolution of singularities has stimulated mathematical research 
for more than a century now, ranging from first results for curves in the late 19th 
century through to the case of surfaces (e.g. [M], [15]) and 3- folds (e.g. [T8]). 
culminating in the famous proof of Hironaka for the general case in character- 
istic zero in 1964 [12]. But this result did not end the interest in the problem: 
First of all, the case of positive characteristic is still an open problem and an 
active field of research; on the other hand, the fact that Hironaka's proof was 
highly non-constructive, led to the question of the existence of algorithmic ap- 
proaches in characteristic zero and to the development thereof (e.g. [2].|16j.[9] 
etc.). 



In its simplest form the problem of desingularization can be phrased as: 
Given a variety X , determine a non-singular variety X and a proper birational 
morphism tt : X — !• X which leaves X \ Sing(X) unchanged. In general, it 
is a rather complicated task to construct such a morphism tt by an iteration of 
improvements of X through blow-ups in appropriate centers, where the choice of 
center turns out to be the crucial point. In certain special situations, however, 
like binomial varieties or curves and surfaces, significant simplifications in the 
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choice of centers and the computation of the blow-ups or the additional use of 
normalization steps can lead to faster algorithms which sometimes also provide 
an X which is covered by fewer charts than the one constructed by the general 
algorithm. This has led to implementations of the various types of algorithms of 
which at least one for each type is available as a Singular library as of version 
3.2. Being offered a variety of algorithms for very similar tasks is certainly a 
nice feature for the experienced user, but novice users are easily overwhelmed 
by it without hints to guide their choice. Therefore a thorough discussion of 
the structural differences and the resulting pros and cons is the topic of this 
article. As it is not possible to explain all these algorithms in detail in one 
article (this would rather require a book) or at least do each of them justice in 
a short description, we only highlight the main ideas and general structure of 
each type of approach and give the corresponding Singular commands and an 
interpretation of the resulting output in explicit examples. 

For comparing algorithms, computer science usually suggests time and space 
complexity considerations. In this case, however, it has been shown in [T] that 
the complexity of the general resolution algorithm is in the complexity class 
£;m+3 Jqj. ^ 771-dimensional varietjQ, which for a computer scientist basically 
translates to 'it is completely hopeless to calculate anything with this algo- 
rithm'. On the other hand, algebraic geometers know from experience with 
Groebner basis computations or primary decomposition (both in that many 
examples, which they want to compute, are quite far from the maximum com- 
plexity of the algorithm. This suggests to select a set of typical examples from 
practical applications and use these for the comparison, which is our point of 
view in this article. 

After very briefly sketching the common main ideas of all Hironaka style 
approaches and giving a Singular example thereof in the section [2l we con- 
sider the structurally much simpler binomial and toric case in the section [31 
In the subsequent section 21 we then consider specialized algorithms for low- 
dimensional varieties and then conclude the article with a comparison of the 
algorithms in section |S1 

The authors are grateful to Santiago Encinas, Herwig Hauser, Ana Bravo, 
Orlando Villamayor, Santiago Zarzuela, Nobuki Takayama, Mohammed Barakat, 
Hans Schonemann and Gerhard Pfister for fruitful discussions, some concerning 
the mathematical side of this article, some related to software development and 
applicability of the respective Singular libraries. 

^To give the reader an orientation about Grzcgorczyk complexity classes, the class 
contains all linear functions, £^ all polynomial functions and all towers of exponential 
functions . 
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2 Algorithms refining Hironaka's approach in 
the general case 

Hironaka's groundbreaking 1964 article [T^ on resolution of singularities in ar- 
bitrary dimension in characteristic zero contains both constructive and non- 
constructive arguments. Seeking better understanding of this proof, these non- 
constructive steps were then studied in detail and filled with explicit algorithmic 
constructions leading to the first algorithmic proofs of desingularization in ar- 
bitrary dimension in the late 1980s/early 1990s (see e.g. [2] [E]), which were 
subsequently refined further (e.g. [9]). Although these approaches may differ 
significantly in various constructive steps they all follow the general structure 
of the original proof. 

To explain all variants of algorithmic proofs in the fiavour of Hironaka in 
detail would require at least a book, so we concentrate on giving a brief im- 
pression of the common main ideas of all these approaches and the resulting 
computational tasks - in the hope of providing the context for the correspond- 
ing examples in Singular. 

Theorem 1 (Hironaka, 1964) For every algebraic variety over a field of char- 
acteristic zero, a desingularization can be achieved by a finite sequence of blow 
ups in suitable non-singular centers. 

From this formulation, we already see that two different algorithmic tasks 
arise here: the blow up, which is a Grobner basis computation and hence does 
not pose any significant problems, and the choice of suitable centers, which is 
the very heart of the construction. For our context, the gist of the choice of 
center can best be formulated by means of a resolution function whose maximal 
locus yields the upcoming center. 

To define the resolution function, we work with the ideal defining the al- 
gebraic variety, and we proceed by induction on the dimension of the ambient 
space. This resolution function, which needs to be both Zariski- and infinitesi- 
mally upper semi-continuous, has the following general structure: 

{invd;invd-i; ■ . ■) 

where each component invi stands for the contribution by an auxiliary ideal in 
ambient dimension i. Each of these components is then composed of an order of 
the ideal (or in top dimension in some approaches the Hilbert-Samuel function) 
and some counting function which ensures normal crossing with the exceptional 
divisors. The construction of the auxiliary ideals is quite delicate and a key step 
in Hironaka's inductive argument. To illustrate the flavour of this construction 
without going into too many technicalities, let us consider a power series w.r.t. 
a main variable z 

f = z'' + ai{x)z''~^ H +ak{x) C C{z,x}. 
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It has order at most k and this order is attained if and only if the order of each 
coefficient ai{x) is at least i. This can also be phrased as 

ordo{f) = k ^ oTd^{{at\aJ ,...af''^^-))>k\ 

where the ideal generated by these powers of the coefficients now plays the role 
of the first auxiliary ideal. Of course, the choice of the main variable, i.e. the 
hypersurface of maximal contact, is crucial in this construction and additional 
difficulties here arise from the fact that such a hypersurface exists only locallyll 

There are two implementations of Villamayor- style algorithms for the gen- 
eral case, both in Singular: [8] and [10] where the latter one is part of the 
standard Singular distribution and is well suited for use with the utilities in 
the additional library for further computations with resolution data. In this 
article, we will only give a syntax example for this second implementation: 

Example 1 



> LIB "resolve . lib" ; // load the library resolve 

> ring r=0 ,x(l . . 3) ,dp; // ring of char with 

// variables x(l),x(2) and x(3) 

> ideal J=x(l) ~2-x(2)*x(3) '2; // Whitney umbrella 

> list re=resolve ( J) ; // compute the desingularization 



> size (re); // the output is a list 
2 // of 2 components 

> re[l]; // list of final charts 

// every chart is given by 
// its defining ring 

[1] : // there is only one final chart 

// characteristic : 

/ / number of vars : 3 

// block 1 : ordering dp 

// : names x(2) x(3) y(l) 

// block 2 : ordering C 

> re [2] ; // list of all charts 
[1] : 

// characteristic : 

/ / number of vars : 3 

// block 1 : ordering dp 

// : names x(l) x(2) x(3) 

// block 2 : ordering C 

[2] : 

// characteristic : 



■^For an accessible and concise treatment of an approach of this flavor and the precise 
definitions of the technical terms used here, we refer the reader to [5]. 



4 



// number of vars 
// block 1 

// 

// block 2 



ordering dp 
names x(2) 
ordering C 



x(3) yd) 



The output is a list consisting of two lists of rings. The first component 
gives the final charts, and the second component is the list of all charts. To see 
the information which is inside a concrete ring or chart it is enough to change 
to the respective ring and proceed as follows: 



// define a new ring 

// set this ring 

// show the information 

// the whole space 
// transform ideal 

// list of exceptional divisors 



> def RR=re [2] [2] ; 

> setring RR; 

> showBO(BO) ; 
==== Ambient Space 
_[1]=0 

==== Ideal of Variety: 
_[l]=y(l)-2-x(2) 
==== Exceptional Divisors: 
[1] : 

_[l]=x(3) 

==== Images of variables of original ring: // morphism defining 
_[l]=x(3)*y(l) // the blowing up 

_[2]=x(2) 
_[3]=x(3) 

In addition, using the additional library it is possible to obtain a picture 
illustrating the hierarchy of charts. We include here a Singular syntax example, 
but omit the picture, as it only contains two charts in this example. (Example 
2 in the next section leads to a more meaningful picture included as figure 1.) 

> LIB "reszeta. lib" ; // library needed for visualization 

> LIB "resgraph. lib" ; // library needed for visualization 

> list idenO=collectDiv(re) ; // collect all the exceptional divisors 

> ResTree (re , idenO [1] ) ; // creates a jpg file with the resolution tree 



As this general approach is obviously quite time and space consuming, it is 
worthwhile to use specialized algorithms wherever those are available. 

3 Combinatorial algorithms for the binomial case 

In the setting of binomial ideals or toric ideals, the situation is much simpler 
than in the general case. This allows a restriction of the methods to combi- 
natorial constructions and arguments. Apart from the classical algorithm of 
star subdivision as toric resolution of singularities (see e.g. standard textbooks 
about toric geometry like |12j). there are combinatorial variants of both main 
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algorithmic approaches to Hironaka resolution, i.e. one in the flavour of Bier- 
stone and Milman [3] and one in the philosophy of Villamayor [4] , [5] . 

Here we will only consider the latter one, which is available as [7] in the 
standard distribution of Singular. It is mainly a combinatorial game on the 
exponents of the generators of the ideal. Roughly speaking, the algorithm looks 
for a smallest order term in the binomials and chooses one of the variables in- 
volved in that term to make an induction on the dimension of the ambient space. 
Then, it looks again for the next smallest order. 

This combinatorial game provides an ideal generated by hyperbolic equations 
that is locally monomial. This ideal needs to be rewritten as a monomial one 
to achieve a log-resolution, see m and for detailsi Another benefit of the 
purely combinatorial approach is the applicability of the algorithm in positive 
characteristic (in contrast to the general algorithms). In a bit more detail, the 
philosophy of the computations can be illustrated in the scheme below: 



Arbitrary char Char Computations 



Input: Binomial ideal — > coefficients, exponents 

type integer exponents=rational numbers 

^ type number 

Output: List of rings •<— coefficients, exponents 
(affine charts) type integer 



Example 2 




> LIB "resbinomial.lib" ; 


// 


> ring r = 2, (x(l. .3)) ,dp; 


// 




// 


> ideal J=x(l) ~2-x(2) ~2*x(3) ~2; 




> list B=binResolution( J) ; 


// 


> size (B) ; 


// 


4 




> B[l] ; 


// 


> B[2] ; 


// 


> B[3] ; 


// 


> B [4] ; 


// 



load the library resbinomial 
ring of char 2 with 
variables x(l),x(2) and x(3) 

compute the binomial resolution 

the output is a list of 4 components 

list of final charts in char 
list of all charts in char 
list of final charts in char 2 
list of all charts in char 2 



The first two components of the output are needed for internal computations 
and also for visualization of charts, since the input for the additional library 
has to be in characteristic zero. Then, for visualization of the resolution tree, 
select the first two components. 

> list L2=B[1] ,B[2] ; 

■'This part of the algorithm is not yet implemented (June 2011). 
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> list idenO=collectDiv(L2) ; 

> ResTree(L2,idenO[l] ) ; 




chart 2 




chart 3 




chart 4 


E:l 




E:l 




E:l 




Figure 1: Resolution tree of Example 2. d denotes the dimension of the center 
of the respective blow up. At each chart, E lists the set of exceptional divisors 
visible in the chart. 



4 Algorithmic resolution in low dimensions 

The other main setting, in which specialized algorithms exist, is the situation of 
low dimensional varieties, i.e. curves and surfaces. As before, these situations 
are significantly simpler than the general case and the first approaches to these 
problems historically precede Hironaka's proof by decades. 
As the singular loci of curves can only be finite sets of points, no special pre- 
cautions to choose appropriate centers need to be taken and a resolution can 
be achieved by simply using the points of the singular locus as the centers of 
the blowing ups. For surfaces, the situation already becomes significantly more 
complicated: the singular loci may be 1-dimensional and may themselves be 
singular. In contrast to Hironaka's approach, the original proof of resolution 
of surface singularities, of which the ideas date back to Jung's 1908 article 
([T3],[T7]), does not only involve blow-ups, but also a different type of birational 
map: normalization. 
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4.1 Resolution of surfaces by Jung's approach 



The general philosophy of Jung's approach is to consider a projective surface 
as a ramified covering of the projective plane. The branch locus, i.e. the dis- 
criminant curve, of this covering can then be desingularized until it only has the 
simplest possible singularities - that is until it is normal crossing and all com- 
ponents are non-singular. Modifying the ramified covering in this way and then 
passing to the normalization of the surface, the only remaining singularities of 
the resulting normal surface lie over the singularities of the resolved discriminant 
curve and are toroidal. These singularities can then be resolved combinatorially. 

In comparison to the projection, the resolution of the branch locus and the 
eventual combinatorial resolution of the toroidal singularities, the most expen- 
sive step here is the normalization. 



Example 3 

> LIB "resjung.lib" ; 

> ring R=0, (x,y,z) ,dp; 

> ideal J=x2+y3; 

> list M=jungresolve ( J , 0) ; 

> M; 
[1] : 

// characteristic : 
// number of vars : 4 
// block 1 

// 

// block 2 



// load the library resjung 

// ring of char with 

// variables x,y and z 

// ideal of the surface 

// compute the resolution 

// the output is a list of rings 

/ / here happens to be only one chart 



ordering dp 
names T(l) 
ordering C 



x y z 



4.2 Beyond the geometric case: 

Lipman's construction for two dimensional schemes 

Like Jung's approach, Lipman's approach makes use of normalization. But here 
the surface is not prepared before normalizing it, instead normalization steps 
may occur at many different points in the algorithm. More precisely, given an 
excellent, noetherian, reduced 2-dimcnsional scheme, this algorithm proceeds by 
first normalizing, then blowing up in the singular locus and then continuing in 
this way by normalizing whenever the resulting scheme at an intermediate step 
is not normal and by blowing up the result of this normalization in its singular 
locus again. 

As normalization is a more expensive computation than a blow up, this 
can potentially slow down the algorithm significantly in comparison to Jung's 
algorithm. On the other hand, this approach provides the advantage of being 



8 



applicable for a wider range than the other one even opening up the possibilities 
for considering arithmetic surfaces. 



Example 4 



> ideal J=x(2) ~2*x(3) ~2-x(l) ~3; 

> list rl=lipmanresolve ( J) ; 

> size (rl) ; 
3 



> LIB "reslipman.lib" ; 

> ring r=0, (x(l. .3)) ,dp; 



// load the library reslipman 
// ring of char with the 
// variables x(l) ,x(2) ,x(3) 
// ideal of the surface 
// compute the resolution 
// there are 3 charts 



If we had done the calculation by hand, wc would have used one normaliza- 
tion and one subsequent blowing up whose center is the only singular point of 
the normal surface. 

As before, wc may pass to any of the charts and have a look at the data there: 



5 Comparisons and Timings 

The following table shows several examples computed with resbinomial, resolve, 
resjung and reslipman, over fields of characteristic zero. The computation times 
t are measured in seconds, and the size means the size of the output in terms of 
number of final affine charts. The subscripts b, s, j and I indicates resbinomial, 
resolve, resjung and reslipman respectively. 

Test results with hardware: Intel core 13 processor 330M (2.13 GHz, 1066 
MHz FSB), 4 GBytc memory. 

The computed test examples show the stability of resbinomial with respect 
to the equations generating the ideal. Since the algorithm is basically a com- 
binatorial game on the exponents, similar equations provide similar size of the 
output- increasing with growing size of exponents. Note that the number of 
charts changes considerably also when passing to higher dimensional varieties. 

In the case of resolve we can see the infiuence of the geometry of the problem. 
It is more efficient than resbinomial if at the beginning is possible to consider a 
permissible center which is bigger than the center selected by resbinomial. 



> def S=rl[l] ; 

> setring S ; 

> slocus ( J_new) ; 
_[1]=1 



// first chart 

// set this ring 

// compute the singular locus 

// non singular 
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Table 1: Test examples 



surface or variety 


tb 


is 




ti 


sizeb 


sizes 


sizcj 


sizei 










1 





12 


1 


3 


1 


X^ 





1 


1 





7 


4 


2 


1 


/y.2 1 T^<3 ryi3 1 ^2 ™5 




1 


1 


1 




14 


5 


3 


X]^3^2'-'^3 '^4 


9 


> 7200 






240 








™2™2 /yi2 1 3 

x^ X2X3, X4 X2 


1 


2 


> 7200 


1 


43 


10 




1 




313 


> 7200 


> 7200 


3 


1863 
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i= time in seconds, size^numbcr of final charts. 



With respect to resjmig, the number of charts remains almost constant for 
the examples of hypersurfaces in ambient dimension three. 

In Lipman's approach, sizei = 1 means that we finish after the first nor- 
malization. This approach is the most efficient in the computed examples in 
ambient dimension 3, but this is mainly due to the rather simple structure of 
the chosen examples, which are very accessible to Grauert's normalization algo- 
rithm, the bottleneck of Lipman's approach. 
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